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ABSTRACT 


A method for sensing retrieval of Aerosol Optical Depth (AOD) was investigated by 
Vincent (2006). This technique is known as the Shadow Method. Using high-resolution 
commercial satellite imagery, Vincent was able to calculate AOD values by measuring 
the radiance of a scene in and out of a shadow. Over the last five years, several 
advancements have been made to validate the Shadow Method. Using the MODerate 
resolution Imaging Spectrometer (MODIS), scenes were analyzed over desert regions to 
exploit shadows generated by clouds and terrain. The results were quickly compiled 
using MATLAB. Results confirm that the Shadow Method 1s capable of producing AOD 


values. 
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I. INTRODUCTION 


The earth’s atmosphere 1s comprised of matter in three physical states: liquid, gas 
and solid. Atmospheric matter ranging from a few nanometers to 100 micrometers in 
diameter, is classified as an aerosol (Waghelstein 2010). Concentrations of aerosols can 


vary greatly with respect to time and the environment, but can reach numerical densities 


as high as 10’—-10° cm” (Seinfeld 2006). These relatively small particles can originate 
from either natural or man-made sources. Naturally occurring aerosols most often appear 
in the form of volcanic ash, dust, smoke, pollen and sea spray. Man-made aerosols are 
most commonly found to include combustion-generated byproducts such as exhaust from 


automobile and power generation facilities. 


With such great variation in concentration and size, it is important to understand 
the impact aerosols have on the environment. Aerosols alter the radiation budget of the 
earth’s atmosphere by scattering or absorbing direct solar radiation. How much they can 
absorb and scatter varies between aerosol types and is most dependent on shape and size 
(Seinfeld 2006). Greater aerosol concentrations can also increase the amount of cloud 
condensation nuclei in a region. This produces more cloud cover, higher cloud 


reflectance, lower visibility and colder atmospheric and surface temperatures 


Unfortunately, aerosol measurement is not an easy task. Dedicated and expensive 
equipment is required to measure aerosol optical depth (AOD), but is impractical for use 
outside of academia and research. With large variations in aerosol content from place to 
place and time to time, these sensors are limited in delivering relevant information and 
usually have small temporal and spatial windows. Also, most countries do not have the 
funds or the interest to invest in dedicated sensors, so obtaining measurements in 


military-sensitive areas 1s a difficult challenge. 


One method that is gaining momentum is Vincent’s Shadow Method (2006). 
Using high-resolution commercial satellite imagery, Vincent was able to calculate AOD 
values by measuring the radiance of a scene in and out of a shadow. Research in the last 


five years has further validated and refined this process. Evans (2007) analyzed the 
| 


output of the technique for various ground surfaces using QuickBird multi-spectral and 
panchromatic sensor channels. Dombrock (2007) developed a method to automate many 
of the calculations the Shadow Method relied on to obtain AOD values more quickly. 
Sweat (2008) showed that shadows from cloud cover could be used to calculate radiance 
values needed for AOD retrieval. Rivenbark (2009) studied the low-bias trend the 
Shadow Method was producing when compared to in-situ observations and determined 
better selection criteria for shadow generators to help reduce the disparity. Belson (2010) 
wrote a sophisticated code that automated the entire process, from user-defined shadow 
selection to AOD retrieval. Waghelstein (2010) departed from the high spatial resolution 
sensors and analyzed scenes taken with the TACSAT-3 satellite, launched on May 19, 
2009, with the onboard Advanced Responsive Tactically Effective Military Imaging 
Spectrometer (ARTEMIS). This hyper-spectral imagery provided a unique opportunity in 


that images were low in spatial resolution, but were high in spectral resolution. 


This research originally started using hyper-spectral imagery, but the Department 
of Defense mandated that classified data could not be transferred to removable media. It 
was then where the decision was made to use Moderate Resolution Imaging Spectrometer 
(MODIS) imagery. According to NASA, MODIS 1s a key instrument aboard the Terra 
and Aqua satellites. Terra's orbit around the Earth 1s timed so that it passes from north to 
south across the equator in the morning, while Aqua passes south to north over the 
equator in the afternoon. Terra MODIS and Aqua MODIS are viewing the entire Earth's 


surface every | to 2 days, acquiring data in 36 spectral bands, or groups of wavelengths. 


I. BACKGROUND 


A. PREVIOUS RESEARCH 


1. Contrast Reduction Method 


One of the first attempts to identify aerosol optical depth by remote sensing was 
done by Odell and Weinman (1975) via the Contrast Reduction Method. This method 
produced theoretical AOD values by using a contrast modulation function, radiance 
values from two distinct objects with different albedo, mean surface reflectance, sun- 
sensor geometry and a specific aerosol model. Kaufman and Joseph (1982) continued this 
technique using a “Two Halves” approach that relied on imagery with obvious 
discontinuity in albedo, like a coastline or forest boundary. Both of these methods were 


heavily dependent. 
pa Dark Object Method 


The Dark Object Method uses surface features with low reflectance to find optical 
depth contributions from aerosols. Any measured radiance beyond the assumed surface 
albedo of a low surface reflectance feature, like a deep body of water or dense forest 
region, 1S attributed to aerosol reflection along the path. Kaufman and Sendra (1988) 
created an algorithm to apply the technique over dense vegetation. Much like the Contrast 
Reduction Method, the Dark Object Method has the tendency to produce large errors due 
to false assumptions in aerosol models and surface albedo. This method also cannot be 


used in areas of high surface reflectance like deserts or urban areas. 


The Dark Object Method was expounded upon when Hsu et al. (2004) developed 
the “Deep Blue” Method, which sought to find AOD values in areas of high surface 
reflectance. Hsu’s team exploited the fact that bright objects with high reflectance often 
reflected back strongly in longer wavelengths, but were not as reflective in the shorter 
wavelengths that corresponded to the satellite’s blue channel. By using sensors with two 
or more blue channels, AOD values could be derived, depending on if the surface 


3 


reflectance was not too high. In the best cases, the Deep Blue Method achieved values 
within 20% of NASA’s Aeronet values, but was limited to platforms with multiple blue 


channels. 
3: Multi-angle Method 


Another approach aimed at finding AOD was the Multi-angle Method, developed 
by Veefkind in 1998. Veefkind conceived a way to retrieve AOD in areas of high surface 
reflectance without the need for multi-blue satellite channels. This method used satellite 
platforms capable of multiple viewing angles collected simultaneously. By collecting 
images at different angles simultaneously, atmospheric columns from each image could 
be categorized and used to derive AOD. While the process proved to be accurate, both in 
high and low surface reflectance environments, it had two major limitations. First, the 
method was restricted to multi-angle sensors. Secondly, these sensors were restricted to 
three- to nine-day revisit times due to their orbital parameters (Veefkind et al., 1998), 


making tactical use over a sustained region of interest unfeasible. 


B. SHADOW METHOD 


1. Introduction 


Unless otherwise cited, the information provided in this description of the Shadow 
Method was taken from Vincent (2006). The Shadow Method exploits the shadows cast 
by buildings, clouds and terrain in high-resolution imagery to measure total optical depth. 
This research sought to exploit shadows cast by clouds and mountains over desert terrain. 
Several commercially-available satellites have sufficient shadow detection capability. 
This research used Moderate Resolution Imaging Spectrometer (MODIS) imagery aboard 
the Terra Satellite. 


In order to use the Shadow Method, it is important to understand the components 
used to get the total and aerosol optical depths. Figure 1 depicts the basic components of 
radiance received at a satellite, both in-shadow and out-of-shadow. Direct transmission 


is the amount of radiation that travels from the sun to the ground, then reflected back to 


4 


the sensor. Diffuse transmission is the amount of radiation that has reached the sensor 
after it has been scattered by some atmospheric component and reflected back to the 
sensor. Diffuse reflection is similar to diffuse transmission, except the radiation goes 
through another scatter by some atmospheric component before it is reflected back to the 
sensor. All three components make up the out-of-shadow signal, while the in-shadow 
signal lacks the direct transmission component. Taking the difference between the shaded 
and unshaded photon budgets, we get the difference in observed radiance as a function of 
surface reflectance, mean aerosol reflectance, solar zenith angle, satellite zenith angle, 
solar irradiance and total optical depth (TOD). The TOD portion of the signal is the key 
to calculating optical depth. 





Direct Transmission O= 0 
mee 


Ll 
ire 

Diffuse Transmission 

Mylot (14, ) 


Diffuse Reflection 
Back to Surface 



















Wal 


Figure 1. Depiction of the difference between the in-shadow and out-of-shadow signal 
received by a satellite. Subtracting the two signals isolates the direct transmission 
component of the signal which can be used to calculate total optical depth (from 
Vincent 2006). 


2. General Calculation 


In order to derive the AOD using the Shadow Method, a governing equation must 
be used to account for the flux densities that contribute to the top of the atmosphere 
(TOA) radiance. Vincent (2206) provides extensive proof in his work, which will be 
summarized as follows: 

The three flux densities in Figure | are: 

Direct Transmission is given as, 

—O 0/40 


oF ve (1) 


where /4o1s the cosine of the solar zenith angle, fois the spectral solar radiant flux 


density and 00 is the optical depth of the atmosphere. 
Diffuse Transmission 1s given as, 
[oF ot (Oo, LL) (2) 
where ¢ 1s the transmittance. 


Diffuse Reflection is given as, 
aLst (3) 


where L: is the surface radiance and r is the mean aerosol reflectance. 


Given these three quantities, the total flux density can be expressed as, 


OO0/L 


= 0 = 
Ls = r( [oF ve + oF ot(Oo, lo) + TLV) (4) 


where r. 1s the surface reflectance. 


Solving for the Surface Radiance, Lx, yields, 


i _ rs Moko 





—0 0/ Lo 
(e + t(So, fo)). (5) 


l-rnr 
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For the satellite to sense the measured surface radiance of an unshaded region of 
the scene, the term must be multiplied by another attenuation coefficient (atmospheric 
optical depth) as it travels through the atmosphere a second time to reach the sensor. 


Thus, the unshaded surface radiance can be expressed as, 


: ofo . Ooi zto —O 0/LL0 
[itd meee tise (7 + (So, Mw) : (6) 
l-rr 


The additional extinction term is dependent on the viewing angle, where yw is the 


cosine of the sensor zenith angle. 


To obtain the measured radiance of a shadow region, the equation, 
0 0/0 


| is = fer Lo Fo ((340))}¢ iH (7) 


l-rnr 4 
is basically the same as Equation (6), with the exception of the direct transmission term. 


Assuming that the shaded and unshaded regions are both a homogenous surface, 
the difference in radiance 1s, 
64 Lt 
Vs LoF 0 ae ad 


L. unshaded — L. shaded = fe = = e : (8) 
l-rnr 








So, solving Equation (8) for Total Optical Depth (TOD) yields, 


5, =| Ho Jin] | —S— ]] Hote |). (9) 
: LU + Ly l=7f htt. 


where 6, is total optical depth, 1, 1s the cosine of the solar zenith angle, is the cosine of 


the sensor zenith angle, r, is the surface reflectance, ris the mean aerosol reflectance 
(MAR), Fo is the spectral solar irradiance for the specific channel, and Lg is the difference 
between the nonshadow and shadow radiance. The equation for calculating TOD 1s 
derived from Liou’s (2002) development of the principles-of-invariance (POI) method to 


describe the multi-scattering problem. 


3. Molecular Rayleigh Scattering Optical Depth Calculation 


In order to derive an accurate AOD measurement from the TOD, a correction is 
applied that accounts for the optical depth due to Molecular Rayleigh Scattering. Russell 
et al. (1993) provides a well documented approach to calculating Rayleigh Optical Depth 
(ROD) values dependent on wavelength, pressure and height: 


Sp (2) = (0.00864 + 6.5x10°H) 4” (p/p,) (10) 
b(A) = 3.916 +0.0742 4+. 0.050/2 (11) 


Where 4H is the height above sea level of the radiometer (km), p is the atmospheric 
pressure at the height of the radiometer (hPa), A is the wavelength (micrometers) and po 1s 


the sea level reference pressure, one standard atmosphere (1013.24 hPa). 
4. Calculating Aerosol Optical Depth (AOD) 


Equation (9) provides a useful tool for measuring TOD. Most of the variables are 
either known and the radiance difference between the shadow and nonshadow region (Lg) 
and the surface reflectance (r,) are measured from the image. Solar and sensor zenith 
angles can be found using the Terascan software suite. The solar irradiance (Fo) 1s 


constant for the specific channel being used, 1621 and 924 Wm’w for Band 1 and Band 


2, respectively. 


Deriving AOD is an iterative approach well outlined in Vincent’s research. 
Basically, once the needed variables are collected, a first iteration of TOD is calculated 
under the assumption the mean surface reflectance (r,) is much larger than the aerosol 
reflectance (r ). This allows the top-of-the-atmosphere (TOA) reflectance to be used as 


the surface reflectance in the shadow region. 


Once the first iteration of TOD 1s obtained, it can be used in the following 


equation to find the Mean Aerosol Reflectance (MAR), 





1 ~ yo = Ene 
(i ES l—e i ‘i du'dd'du"d¢" (12) 
0 


where @ is the single scatter albedo, dis the TOD and P(®) is the scattering phase 


function. 


The scattering phase function in this case is the Henyey-Greenstein phase 

function, 
p(6)=_—__-2 _. (13) 
(1+ 2° —2gcos®) 

where © is the scattering angle and g is the asymmetry parameter that can range from -1 
to 1, fully backscattered condition to fully forward scattered, respectively. Values fro the 
single scatter albedo and the asymmetry parameter can be assumed based on expected 
aerosol type, although Vincent shows that the AOD values do not vary greatly with 
changes to @ or g. Total uncertainty for the Shadow Method is typically +0.04 optical 
depth units based on the assumption in @ and g, but in extreme cases can be as high as 
0.10 optical depth units. The MAR is determined using a numerical approach to the 
quadruple integral in Equation (12). MAR approximates the diffuse sky radiance. 


Once the MAR is calculated, it can be subtracted from the top of the atmosphere 


reflectance to solve for the mean surface reflectance. Since r,and r are known, a second 
and final iteration of TOD can be calculated. Once the final TOD 1s calculated, the ROD 
(Equation 10) is subtracted and we finally get a value for AOD. 


C. MODIS 


Previous research by Waghelstein used hyper-spectral imagery from TACSAT. 
There were several limitations. The first was that the data was classified. The second was 
that most of the satellite passes were over areas nowhere near NASA Aeronet sites. Three 
locations happened to be within 10km, but were urban areas. While urban areas provide 


good shadows, the problem is that the surfaces were not all homogeneous. Shadows cast 
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by buildings encompassed other buildings, roads, houses, etc. This research focused on 
desert areas that had shadows cast by clouds and terrain. The cases used definitely had 
homogeneous surfaces, but one drawback to using desert regions is that other than a dust 
storm, the aerosol content is low. Initial research was focused on finding HS imagery, but 
classified data could not be transferred to removable media and ordering images required 


Artemis to mail data discs. It was then where the focus switched to using MODIS 


imagery. 
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Hl, METHODOLOGY 


A. OVERVIEW 


The Shadow Method was developed with the intent of exploiting regions of 
shaded and unshaded radiance found in high resolution imagery. Vincent stated that the 
image would have to be orthorectified and calibrated for the Shadow Method to retrieve 
accurate AOD values. Follow-on research by Rivenbark (2009) would later prove that the 
time-consuming orthorectification process played no part in the constant low-bias and 
could therefore be ignored in pre-processing stages. Shadow and non-shadow pairs were 
manually collected by the user in the Environment for Visualizing Images (ENVI) 
version 4.7 software suite and run through a customized code in Mathwork’s Matrix 


Laboratory (MATLAB) software that calibrated the image and calculated AOD values. 


B. SATELLITE IMAGERY 


1. Imagery Retrieval 


Each case in this research was obtained from the United States Geological Survey 
(USGS) through the NASA Land Process Distributed Active Archive Center (LP 
DAAC). The files can be ordered through their website at https://Ipdaac.usgs.gov. 


Zs ROI Sampling 


Using the ENVI software suite, Regions of Interest (ROI) were sampled. The user 
would load the appropriate image file and identify regions of shadow/non-shadow pairs. 
It is very important that the user choose to analyze the image in radiance, not reflectance. 
The ROI statistics function allows the user to isolate regions in the scene. Once the region 
is selected, the user then views the “stats” function. There, the user can see the maximum, 
minimum, mean and standard deviation values of brightness values for each pixel for the 
particular region. The user must save the data as a text file in order for MATLAB to read 


it. As proven by Dombrock (2007) and Evans (2008), AOD values were most accurate 
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when compared to the mean calibrated radiance value within each ROI pixel, as opposed 


to a minimum, maximum or mode sampling values. 


In order to pick the shadow/non-shadow pairs that would yield the best results, 
care must be taken by the user. Dombrock (2007) and Rivenbark (2009) showed that 
shadow ROIs selected within the center of the shadow produced better results than 


selecting large portions of the shadow or near the edge of the region. 
aP AOD Calculation 


The MATLAB code was originally written in the Interactive Data Language 
(IDL) for ENVI by Vincent (2006) and later altered by Dombrock (2007) and Rivenbark 
(2009). In 2010, Waghelstein translated the code into MATLAB language due to the 
spectrally robust ARTEMIS Hyper-Spectral sensor and the software’s matrix math 
capability. Since MODIS has only two bands compared to ARTEMIS’s 480 bands, the 
code was manipulated once again. Waghelstein’s entire code took roughly 50 minutes to 
compile results. Using the MODIS sensor’s two bands, the code only takes about three 


seconds. 


In order to “simplify” the code, multiple changes were made. Due to the 480 
bands that Wagelstein used, each wavelength for each respective band was pulled from 
the Whereli excel file. Using MODIS, the central effective wavelength could just be set at 
0.650 and 0.860 for Band | and Band 2, respectively. The scale factor, a variable for the 
aperture radiance calibration, was made a constant because the images were already in 


terms of radiance. 


It is recommended that the user identify as many pairs as possible. With multiple 
ROIs selected for each image file, the user is able to save the statistics output into an 
ASCII text file. Using the “total shadowtype.m’” MATLAB file, the ASCII text files are 
imported into the code that would calculate AOD values for each shadow/non-shadow 


pair. 


Since editing the code can get cumbersome, it is recommended the user make 


separate .m files to isolate the type of shadow that 1s being derived, for example, cloud 
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and terrain. Before the code can be run, the user is required to find the solar and zenith 
angles using the Terascan software suite, the number of pairs being investigated, the 
Julian day and ground truth for both bands. Once the code is run, results are displayed in 
nearly three seconds in terms of the AOD for each shadow/non-shadow pair, separated by 


wavelength band. 
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IV. RESULTS 


A. | FT. HUACHUCA, DECEMBER 23, 2009 


wi 
- 
4 


: 
’ 


| Hh ) j 
aay, 


bs 
rf 


t i 
Ba 
; 


& 
ow 





Figure 2. MODIS Image of Fort Huachuca, December 23, 2009. All terrain and cloud 
pairs were taken from the annotated area. 
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1. Overview 


The first case analyzed was at Fort Huachuca on December 23, 2009. This site 
was originally chosen using the Hyperspectral imagery from the ARTEMIS sensor 
because of its close proximity toa NASA AERONET sensor. The image was ordered, but 
could not be mailed due to classified restrictions. The same image could have been 
ordered from the MODIS site. This scene was ideal because it had both cloud and 
topography shadows. All of the pairs were taken from the northern section of the image, 


along the cloud border (Figure 2). 
2 Hypothesis 


The original hypothesis was that all Shadow Method-derived AOD values, both 
cloud and terrain-generated, would be within 0.05 units of the NASA AERONET ground 
truth data. I also expected the range of values to be small, within 0.02 units. Since the 
entire scene is fairly homogenous, cloud and terrain-generated shadows were used that 
had shadows larger than 30 pixels. Shadows any smaller than 30 pixels would be difficult 
to analyze using the ROI tool. 


3: Cloud 


In all result figures, each symbol represents the Average AOD for a particular 
shadow/non-shadow pair. The values for the cloud pairs range from close to zero to 0.06 
for Band | and 0.03 to 0.09 for Band 2 (Figure 4). Band 2 is higher in all cases due to the 
mathematics. The Central Effective Wavelength (CEW) for Band | and Band 2 1s 0.650 
and 0.860 respectively. This makes the Rayleigh Optical Depth (ROD) larger for Band 2. 
The solar irradiance (Fo) is 1621 and 924 Wm’ for Band 1 and Band 2, respectively. 


(Fo) is in the denominator of the “albedol” function in the code, making Band 2 larger. 
The ground truth from the AERONET sensor was 0.095 and 0.10 for Band 1 and 
Band 2, respectively. These are very low values for AOD. This was expected due to the 


location of Fort Huachuca, which is far from any urban area and apart from clouds, was 


clear with respect to smoke, dust and pollution. 
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Pairs 2, 4 and 6 yielded the results closest to the ground truth. This is due to 
greater contrast in radiance between the shadow and non-shadow area. These pairs were 
within the 0.05 AOD unit expectation. Pair 1 appeared to have white sand in the scene 
(Figure 3A). This area of higher albedo has a higher reflectance, lessening the contrast 
between the shadow and non-shadow region. Pair 3 incorporated a very small shadow of 
24 pixels (Figure 3B) and pair 5 was surrounded by cloud (Figure 3C). Higher resolution 


imagery could have produced more accurate results. 


Pair 7 1s consistent with what Dombrock (2007) and Rivenbark (2009) showed in 
their research where ROIs selected within the center of the shadow provided better results 
than selecting large portions of shadow near the base of the generator or the edge of the 
shadow/non-shadow region. The shadow region selected in this scene was very close to 
the border (Figure 3E). The differences are rather small, but considering the AOD was 
around 0.10, values above zero using the shadow method are encouraging. The level of 
uncertainty with the range of values was perplexing, even though it is consistent with 
what Vincent showed with total uncertainty for the Shadow Method is typically +0.04 
optical depth units. Ideally, the AOD values should have been analogous. The regions 
that were sampled were taken in close proximity to one another and the user assumed that 
the atmosphere was homogenous. It was concluded that the variation in the results were 
caused by the shadow generator or the characteristics of the surface. The scenes that were 
selected appeared to be mostly homogenous, with only minor differences in the terrain. 


This led to the analysis of regions with terrain-induced shadows. 
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Figure 3. Cloud pairs from Fort Huachuca, December 23, 2009. A) Pair 1 (top) and Pair 
2. B) Pair 3 (top) and Pair 4. C) Pair 5. D): Pair 6. E) Pair 7. 
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Cloud: Terra, 12-25-2009, 2055 UTC, 1225 
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Figure 4. MATLAB results from Fort Huachuca Cloud Pairs, December 23, 2009. Pair 
: red square, Pair 2: blue square, Pair 3: green square, Pair 4: aquamarine square, 
Pair 5: red circle, Pair 6: blue circle, Pair 7: green circle, Pair 8: aquamarine 
circle, Computed Average AOD: black diamond, NASA AERONET AOD 
(Ground Truth): Pink Star. 
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4. Terrain 


The AOD values from the terrain cases were not as close to the ground truth as 
the cloud results (Figure 6). Only the fourth pair was within the 0.05 AOD unit 
expectation. However the range of values is much smaller. This is due to characteristics 
of the shadow generators. The terrain is homogenous, compared to the cloud generators. 
With cloud, radiation can get through, allowing some to reach the ground where it could 
be reflected up to the sensor. With terrain, radiation cannot get through the shadow 
generator, making the shadow a “true” shadow. The low AOD values, compared to the 
cloud results could be due to the size of the shadow pairs. The average number of pixels 


for each pair was about 32, compared to the cloud scenes of 100. This limited the user 
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from avoiding the edge of the shadow. Dombrock (2007) and Rivenbark (2009) showed 
that shadow ROIs selected within the center of the shadow produced better results than 
selecting large portions of the shadow or near the edge of the region. The results were 
much more grouped and Band | was lower, as expected. The values for the results in 
Band 2 were higher, but still less than half of what the AERONET sensor recorded. Pair 4 
yielded the lower results due to the snow near the shadow region (Figure 5C). Just like 
with the cloud case with pairs 1, 3 and 5, higher albedo in the shadow region seem to 
reduce the contrast between the shadow and non-shadow regions. The area of snow was 
analyzed as the non-shadow, but the area was only 20 pixels. The range of values in this 
case was more impressive due to the tight grouping, but not convincing for accurate 


Shadow Method-derived AOD values. The ground truth value of 0.10 essentially shows 


that the air was clear and no significant value of AOD could be obtained. 





D E F 


Figure 5. Terrain Pairs from Fort Huachuca, December 23, 2009. A) Pair 1 (top) and 
Pair 3. B) Pair 2. C) Pair 4. D): Pair 5. E) Pair 6. F) Pair 7. 
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Terrain Terra, 12-24-2009, 2055 UTE, 1225 
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Figure 6. MATLAB results from Fort Huachuca Terrain Pairs, December 23, 2009. Pair 
1: red square, Pair 2: blue square, Pair 3: green square, Pair 4: aquamarine square, 

Pair 5: red circle, Pair 6: blue circle, Pair 7: green circle, Computed Average 

AOD: black diamond, NASA AERONET AOD (Ground Truth): Pink Star. 


21 


B. SAUDI ARABIA, DECEMBER 12, 2003 





Figure 7. MODIS Image of Saudi Arabia, December 12, 2003, 0710 UTC. All cloud 
pairs were taken from the yellow-annotated area. 
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1. Overview 


After seeing the results from the Fort Huachuca case, it was decided to analyze 
cases that were known to have higher AOD values. A dust storm over Saudi Arabia on 
December 12, 2003 was chosen (Figure 7). The image was taken at 0710 UTC, which is 
1140 local time. Shadows are hard to come by at this time because the sun is almost 
directly overhead at noon. Shadows were found, but mostly in the northern portion of the 


image. Shadows generated by terrain could not be found because the region 1s mostly flat. 
2; Hypothesis 


Since ground truth could not be obtained, A different approach had to be taken. 
Since this is a region of high AOD, it was desired to retrieve Shadow Method-derived 
values larger than the previous case. Even though the area that was sampled was in a 
clearing area, the values should be larger than 0.10. The previous case showed that results 
vary on the characteristic and size of the shadow generator. The goal of this case was to 
show that the Shadow Method could produce results larger than 0.10. With respect to the 
range of values, 0.05 was acceptable like the previous case since not all clouds are the 


same. 


3. Cloud 


Clouds were present on the southeastern edge of the image, but were too small to 
try, roughly 10 pixels. Noon is early in the day and the clouds were just beginning to 
form. After looking over the entire image, shadows large enough for using the ROI tool 
were found over Southwestern Iran. Four pairs were analyzed that yielded an AOD 
average of 0.065 (Figure 9). These values seem low and make sense due to the clearing 


area in the north. Similarly to the first case, the results were higher for Band 2. 


Pair | has a lower value due to the thickness of the cloud (Figure 8A). The cloud 
appears to have a lower brightness temperature. This means that the cloud is warmer in 
temperature and not as thick, allowing more radiation to reach the ground. The cloud is 


thick enough to cast a shadow, but the radiance difference is not as great as the other 
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three pairs. Pairs 2, 3 and 4 clearly show brighter clouds, indicating higher cloud tops. 
Since these pairs were taken fairly close to each other, it 1s assumed that the Lifted 
Condensation Level, the level where clouds form, is the same. If the cloud base level is 


the same, then brighter cloud tops would indicate thicker clouds. 


The take-away from this case is encouraging because the Shadow Method does 
yield results, but just like the Fort Huachuca cloud case, they do vary based on the 
characteristics of the shadow generator. With the exception of the first pair, I was pleased 
with the values and range of values. Pairs 2, 3 and 4 were within 0.02 of each other, 
exceeding my expectations. If the image was taken later in the day, perhaps terrain- 
generated shadow scenes could have been analyzed. The more shadow pairs the user can 
find, the more accurate the AOD average will be. Since four pairs were analyzed here, the 
range of results was quite large. Ground Truth could not be obtained due to the fact that 
most countries, like Iran, are not forthcoming with their weather data or cannot afford the 


technology to share their data. 








Figure 8. Cloud Pairs from Southern Iran, December 12, 2003. A) Pair 1. B) Pair 2. C) 
Pair 3. D): Pair 4. 
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Cloud: Terra, 12-12-2003, OF 10 UTC, 34° 
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Figure 9. MATLAB results from Saudi Arabia Cloud Pairs, December 12, 2003. Pair 1: 


red square, Pair 2: blue square, Pair 3: green square, Pair 4: aquamarine square, 
Computed Average AOD: black diamond. 
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C. IRAN/AFGHANISTAN, DECEMBER 11, 2003 


KO 





Figure 10. MODIS Image of Iran and Afghanistan, December 11, 2003, 0625 UTC. All 
cloud pairs were taken from the yellow-annotated areas. Terrain pairs were taken 
from the red-annotated areas. 
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1. Overview 


The third and final case analyzed was the same region as the second case, but 23 
hours later. The satellite pass at this time was also further northeast. This image is ideal 
because it had more shadow generators, both due to cloud and terrain. These areas also 
varied in that they lay in areas of high and low AOD due to the dust storm (Figure 10). 
Shadows could not be analyzed in the thicker areas of dust. As with the previous case, 


ground truth could not be obtained. 
2. Hypothesis 


Since this scene has varying regions of AOD, the goal was to show that areas in 
the vicinity of the dust storm could produce distinctively higher Shadow Method-derived 
AOD values for both the cloud and terrain-generated pairs. Using 0.10 as a minimum 
value made sense, since it was used in the previous two cases. The range of values, 
however, needs to be looked at differently, since the scene varies. A large range of values 
for both the cloud and terrain-generated shadow pairs was acceptable because a large 


portion of the scene was sampled. 
3. Cloud 


This case produced higher results, compared to the previous case. Values ranged 
from 0.07 to 0.29, the highest results of all three cases (Figure 12). Pairs 1-3 and 6—9 
yielded the lowest results. Taking a closer look at Figures 11E-—F, it is obvious that pairs 
6 and 7 are clearer than the others. This is also evident in Figure 10. There is some dust 
over the scene, but not much. These values are somewhat close to the values for pairs 1—3 
and 8—9, even though there 1s less dust. The reason for this is similar to the previous case 
in that cloud characteristics play a role. In Figures 11E—F, the brightness temperature is 
higher than the others. This indicates thicker clouds, which means that the cloud is 
allowing less radiation to reach the surface. So, the shadow is darker for pairs 6 and 7, 
making the difference between the shadow/non-shadow regions greater, giving a higher 
value for AOD. If the clouds were all the same, then pairs 1-3 and 8-9 would show 


greater AOD values, since there is more dust over those scenes. 
2) 


Pairs 4 and 5 yielded the highest results (Figure 12) because this region was 
surrounded by the dust storm (Figure 10). The brightness temperature for this scene 1s 
lower, but the dust content is greater, which won out. If the clouds were as thick as pairs 


6 and 7, the AOD values would have been even greater. 
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Figure 11. Cloud Pairs from Iran and Afghanistan, December 11, 2003. A) Pair 1. 
B) Pair 2. C) Pair 3. D): Pair 4 (left) and Pair 5. E) Pair 6. F) Pair 7. =G) Pair 
8. H) Pair 9. 
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Figure 12. MATLAB results from Iran/Afghanistan Cloud Pairs, December 21, 2003. 


Pair 1: blue square, Pair 2: green square, Pair 3: aquamarine square, Pair 4: 
red circle, Pair 5: blue circle, Pair 6: green circle, Pair 7: aquamarine circle, Pair 
8: red star, Pair 9: blue star, Computed Average AOD: black diamond. 


1. Terrain 


Three pairs of terrain scenes were analyzed in Afghanistan. These areas were not 
in the thick portions of the dust storm (Figure 10). The images appear as clear as the 
cloud pairs in Figure 11E—F. The values for those cloud pairs (6 &7) in Band 2 ranged 
from 0.10 to 0.125 (Figure 13) and the values for this terrain case ranged from 0.05 to 
0.125 in Band 2 (Figure 14). Even though the cloud and terrain pairs were analyzed in 
different areas, they should exhibit similar AOD values. 


Pair 2 exhibited the higher results due to the scene being on the northern edge of 


the dust storm (Figure 10). The area was clearing, but aerosols were still present. Pairs | 
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and 3 are on the eastern edge of the storm and judging by the northwesterly wind 
direction, have not been greatly affected yet. If the image were taken an hour later, AOD 


results would have been higher. 


A B C 


Figure 13. Terrain Pairs from Afghanistan, December 11, 2003. A) Pair 1. 
B) Pair 2. C) Pair 3. 
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Figure 14. MATLAB results from Afghanistan Terrain Pairs, December 21, 2003. 
Pair 1: red square, Pair 2: blue square, Pair 3: green square, Computed 


Average AOD: black diamond. 
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V. CONCLUSION 


A. FINDINGS 


All three cases provided varied results. Fort Huachuca was analyzed on a very 
clear day. The results were low, but greater than zero, which was encouraging. The 
cloud-generated shadow pairs theoretically should have been similar to one another, but 
were varied depending on the characteristic of the cloud. One would not believe any 
worth can come from cloud-generated shadow analysis. Values can be retrieved, but how 
reliable would they be? Not all clouds in a scene are the same, so one would get varied 
AOD values in an area that should get the same. The Afghanistan case is a good example 
of this. In an area of thick dust, the Shadow Method-derived AOD for cloud pairs 8 and 9 
should have much higher, but was low due to the cloud thickness. One would need to 


know the cloud depth and be certain that the shadow is genuine. 


For the Fort Huachuca case, the terrain pairs could have been higher if the shadow 
regions were larger, but all seven pairs had nearly identical AOD values making the 
uncertainty low, but accuracy was also poor. Perhaps if the AOD was higher, more 


convincing results could have been found. 


Both the Saudi Arabia and the Afghanistan cases proved that the shadow method 
does, in fact produce results. Apart from different cloud characteristics, both cases 
yielded results that made sense. Areas of intense dust had higher AOD results than the 
areas that were clear. The terrain case for Afghanistan produced similar AOD values and 
it is supported that both this case and the Fort Huachuca terrain case produced more 


trustworthy results. 
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B. UNCERTAINTY 


While the cloud cases provided the higher results for all three scenes, the results 
were varied based on cloud characteristics. Scenes that should have yielded higher AOD 
values like the Iran and Afghanistan cases ended up producing lower results due to 
decreased cloud thickness. If one were to try to retrieve AOD using the shadow method, 
care must be taken by the user. Areas of high AOD could be underestimated if the clouds 


are not thick enough. 
C, RECOMMENDATIONS AND FOLLOW-ON RESEARCH 


For the next round of research, two avenues are recommended. Since this research 
showed that the terrain-generated shadow scenes produced more consistent results, the 
first recommendation would be to find cases that have high AOD scenes with many 
terrain-generated shadows. Ideally, these scenes should be close to a ground truth source. 
The time of day should also be taken into consideration. Images taken around noon 
would provide smaller shadow regions, so early morning or mid-afternoon would be 


ideal. 


The second round of research would be to look into what cloud depth would 
provide the more accurate AOD. The three cases in this research had more cloud pairs 
than terrain pairs, so a greater number of cases could be analyzed. In order to do this, one 


would need ground truth and upper air data to confirm atmospheric conditions. 
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